Probing Space-time Foam with Photons: Suppression of Observable Effects due to 

Uncertainty in Optical Paths 
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It was recently proposed to use extra-galactic point sources to constrain space-time quantum 
fluctuations in the universe. In these proposals, the fundamental "fuzziness" of distance caused 
by space-time quantum fluctuations have been directly identified with fluctuations in optical paths. 
Phase-front corrugations deduced from these optical-path fluctuations are then applied to light from 
extragalactic point sources, and used to constrain various models of quantum gravity. In partic- 
ular, the so-called random-walk model has been claimed to be ruled out by existing astrophysical 
observations from the Hubble Space Telescope. However, when a photon propagates in three spatial 
dimensions, it does not follow a specific ray, but would rather sample a finite, three-dimensional 
region around that ray — thereby averaging over space-time quantum fluctuations all through that 
region. We use a simple, random-walk type model to demonstrate that, once the appropriate wave 
optics is applied, the averaging of neighboring space-time fluctuations will cause much less distor- 
tions on the phase front. In our model, the extra suppression factor due to diffraction is Ip/X, which 
is at least 19 orders of magnitude for astronomical observations. 

PACS numbers: 95.75.Kz, 03.67.Lx, 04.60.-m, 98.62.Gq 



It was recently proposed to use extra-galactic point 
sources to constrain the quantum fluctuation in space 
time P, 0, EI 0- It was argued that, space-time fluc- 
tuations cause random phase shifts of photon, and that 
these fluctuations accumulate throughout the very long 
light propagation path from the point source to the 
earth, causing wavefront distortion from a perfect spher- 
ical shape upon arrival at the earth. The manner these 
fluctuations accumulate depend on the specific model of 
quantum gravity, and in particular, it has been claimed 
that the random walk model could be ruled out by exist- 
ing imaging data from the Hubble Space Telescope. In 
this paper, we point out what we think is a serious omis- 
sion in the theory so far employed by the above authors, 
and argue that the random walk model, once given a 
closer look, cannot be ruled out at all by current or any 
forseeable observations of extra-galactic point sources. 

In Rcfs. the authors assumed photons origi- 

nating form the point source to undergo a random phase 
shift A(j> due to space-time fluctuations. According to 
various models of quantum gravity, they have 



A<£~ 2tt(; p /A) q (L/A) 1 



(1) 



where lp is the Planck length, A is the wavelength of 
the light, L is the light propagation length, and a is a 
parameter that depends on the model of quantum grav- 
ity. In particular, a = 1/2 corresponds to the so-called 
random-walk model, which can be understood as having 
the speed of light fluctuating dramatically at the scale 
of Zp. The relation Q stems from a more fundamental 
ansatz of quantum gravity, in which the uncertainty in 
distance measurement 8L over a distance L is given by 



which was first derived by Ng and van Dam and then 
discussed by others 0, H, H 13 . Note that Eq. (JTJ is 
related to Eq. (2) by A<j> ~ 2ir(SL/\). 

While Eq. J2J| could be regarded as fundamental to the 
quantum gravity theory, Eq. cannot, because it also 
depends on how light propagates: when a photon trav- 
els through a space-time region, it does not follow only 
one particular ray, whose length is subject to the funda- 
mental "fuzziness" prescribed by Eq. (J2J, but instead, it 
would simultaneously sample an ensemble of many dif- 
ferent neighboring rays, each of which having a poten- 
tially different realization of the fundamental length fluc- 
tuation; the actual path-length fluctuation must then be 
given by an averaging among these different length fluc- 
tuations. This allows A(j> to go below 2irSL/\. More- 
over, because the linear size of the sampling region can 
be much bigger than Zp, the correlation length of funda- 
mental quantum fluctuations, this averaging can dramat- 
ically suppress the actual Acj> from 2ir5L/\, or Eq. 0). 

But we note that, if we lived in a space-time with one 
time dimension plus one spatial dimension, then there 
would be no uncertainty in optical paths, and there would 
not be any suppression factor from 

Model of space-time fluctuation. In the rest of this pa- 
per, wc will elaborate the above argument using a sim- 
ple model of space-time foam. We consider the light as 
a scalar wave traveling in Minkowski space-time, filled 
with a medium with random (yet static) distribution of 
refractive index n(x, y, z) = 1 + e(x, y, z). We suppose 
that e has the following translational invariant spatial 
auto-correlation function: 



SL > L 1 -"^ , 



(2) where "(. 
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unity, and 
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1 \x\ < Ip 
|x| > l P 



(4) 



In this toy model, the coordinate speed of light propaga- 
tion fluctuates strongly, if we zoom in on a small region 
with size comparable to the Planck scale — this simu- 
lates light propagation in a space-time with Planck-scale 
quantum fluctuations. [Note that we have ignored tem- 
poral fluctuations.] In addition, light-speed fluctuations 
in regions separated by more than the Planck length are 
independently from each other. 

For a moment, if we took this mental picture of speed- 
of- light fluctuation, forgot about the wave nature of light, 
and assumed that each photon would follow a distinct co- 
ordinate path in (t, x, y, z), or if we assumed a space-time 
with only one time dimension plus one spatial dimension, 
we could also derive Eq. Q with a = 1/2. In order to do 
so, we divide the distinct coordinate path into Nl = L/lp 
intervals, each with length lp. The phase-shift fluctua- 
tion in each interval is S(f> ~ 2nlp/\, while fluctuations 
in different intervals are independent from each other. 
In this way, the total phase shift of the photon does a 
"random walk" while the light propagates. At the end of 
propagation, we have a photon-phase fluctuation of 



(A0) 



ID 



N L 5cf> ~ yflpL/\* , 



(5) 



which is exactly Eq. (JJJ, with a = 1/2. 

The Wave Equation. Returning to the wave picture, 
we first write down the wave equation: 



<9 2 $ 

[l + 2e(x,y,z)} — + V 2 $ = 0. 



(6) 



Since our refractive-index perturbation is static, we can 
expand the total wave into two monochromatic pieces, 
the primary wave, &o(x, y, z)e~ lu)at and the secondary 
wave, ip( x i z)e~ tuJot . At leading order in e, wc have 

(V 2 + uJq) ip(x, y, z) = -2wq£(x, y, z)$ (x, y, z) . (7) 



For a point source, we assume 



r = y/ ' x 2 + y 2 + z 2 



(8) 



At a distance L, the secondary wave if> must be compared 
with the primary wave at this point to give the amount 
of modulation it induces. We define 



AirLipe' 



-iuj L 



(9) 



with a, 4> G R- If ct 2 + 4> 2 <C 1, a is the relative amplitude 
modulation, and <p is t ne phase modulation in radians. 
We also define the total modulation, 



f = y/(a 2 + cj) 2 ) = 47rLvW*> 



(10) 



which is larger than the standard deviations of both the 
amplitude and the phase modulation. 

Summing Over Paths. To arrive at the answer quickly, 
we use the Huygens-Fresnel-Kirchhoff scalar diffraction 
theory, which is equivalent to applying the outgoing 
Green Function, and obtain M : 
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Note that we have only considered contributions from 
fluctuations at distances smaller than L to the point 
source. The integral l|ll(l can be interpreted as a path 
integral — over all paths that consist of two straight 
sections (each associated with a propagator), and a de- 
flection in the middle due to interaction with refractive- 
index fluctuations (associated with a coupling coeffi- 
cient). Paths with more than one deflection do not have 
to be taken into account in our linear treatment. 

If we discretize the integration domain, a sphere with 
volume ~ L 3 , into cells with linear size ~ lp and volume 
vp ~ lp, we will get a total of iV tot ~ L 3 /lp individual 
cells, each of has a statistically independent fluctuation 
in e with variance a 2 [Cf. Eq. 0]. We will then estimate 
that Eq. i|ll|) would give 

VW) - ^vp ~ aujl ^/L . (12) 

According to Eq. i|l(J|) and comparing with Eq. J3J), we 
have 



>(WA). 



(13) 



There is an extra suppression factor of Ip/X, which arises 
from the fact that in Eq. i|ll|) . the intermediate point 
x' of the optical path has the freedom to go away from 
the axis connecting the source point and the field point, 
and sample through AT to t ~ L 3 /lp independent fluctua- 
tions, instead of only Nl ~ L/lp in the one-dimensional 
treatment. 

A more rigorous calculation using Eqs. J2J, <fTT7f> and 
(|ll|l only gives a more precise numerical factor: 



(14) 



Spatial-Scale Cut-off. To study fluctuations at differ- 
ent spatial scales, we solve the same problem by decom- 
posing the secondary wave into modes: 

*l,(r,6,cf>) = Y^bPem(r)Y em (6,Lp)}. (15) 

Irn 

Here Ye m (0,<p) are Spherical Harmonics. They describe 
angular variations at scales of 2tt/£] at a radius r, that 
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FIG. 1: Plots of £f as functions of £/(ujoL), for cases with 
ljqL — 50 (dotted curve), 100 (dashed curve), and 200 (solid 
curve). All modes with < £/{ojqL) < 1.2 are shown. 



corresponds to transverse length scales of 2wr/£, or trans- 
verse spatial frequency of £/{2irr). The modal decompo- 
sition of Eq. Q is 



1 d 
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(16) 



with 



dip I sm6d6[e(r,6,<p)Y; m (6,<p)] , (17) 



which satisfies 



(ee m (r)eg, m ,(r')) = a 2 5 u >5 mm '6(r - r')l 3 /r 2 
Here we simply assumed 

(e(x)e(x')) = a 2 /jU (3) (x-x'). 



(18) 



(19) 



The spatial spectrum corresponding to this correlation 
function is identical to that in Eq. at low spatial fre- 
quencies, but continues to exist in orders higher than 
1/lp. In principle, those modes may also add incoher- 
ently to our output fluctuations, but as we shall see, their 
contributions will be negligible. Solving Eq. (|16f) . assum- 
ing regularity at r = and outgoing condition at r = L, 
we obtain, 



L>%Lhf~\ujQL) 



2vr 



dr[rj { (u; Q r)e^ r e em (r)}, 
(20) 

where jg and h!jp are Spherical Bessel function and first- 
kind Spherical Hankel function [{J. 
From Eq. {TSJ and we obtain 



167T L {famlptm,') =£ im 6tt>6r, 



(21) 



with 



4a 2 



w Lh { p(L) (loqIp) 3 jj{u r)uj Q dr . (22) 



Note that £f m is independent from m, which is a conse- 
quence of the rotation invariancc of the refractive-index 
fluctuations. The total fluctuation at r = L will then be 



= 16tt 2 L 2 ^(L, 9 9,<p)) 

+oo +1 

= E E ZLYt m (8,<p)Yt m (9 t <p) 



£=0 m= 

-\-QG 



e=o 



1=0 



(23) 



Physically, £ 2 describes fluctuations at the angular scale 
of ~ 2tt/£, or transverse spatial scales of 2irL/£, or trans- 
verse spatial frequency of £/(2ttL). Inserting Eq. i|22[l . 
we have 



a 2 (LO l P f 



2£+l 



uoLh^ (uj L) 



uj L 



]j{R)dR (24) 



We expect £ ~ wo-t, or 2irL/£ ~ A to be the turning 
point, because at this point the transverse spatial scale 
is comparable to the wavelength A. 

Mathematically, for £ < u>oL. the Spherical Bessel and 
Hankel functions are wavelike at r ~ L, indicating prop- 
agating waves; for £ > ujqL, the Spherical Bessel and 
Hankel functions are not wavelike at r ~ L, indicating 
evanescent waves. In the limiting regimes of £ <C ujqL 
and £ 3> <jJ L, ^ can be evaluated analytically, using 
asymptotic expansions of Spherical Bessel functions: 



,c2 



a 2 (uj lpY 



1/2. 



£ < lo L , 



lo L/[(2£+1) 2 tt], £^>u L. 



(25) 



Note that not only does £f approach at orders £ 
luqL, the summation of all these higher modes also gives 
a negligible contribution. This qualitatively confirms a 
cut-off at the transverse scale of A: fluctuations at much 
finer scales do not generate secondary wave. This justifies 
our original use of Eq. I|19|) ■ 

In Fig. ^ we study the transition zone of £ ~ ujqL nu- 
merically, for moderately large values of ljqL = 50, 100 
and 200, by plotting £| as a function of 1/(ojqL). As 
loqL — > +oo, £f asymptotes to a smooth, non-zero func- 
tion for £/{uj L) < 1, and to for £/{lu L) > 1. This 
means, in the realistic situation of u>qL <C 1, £/(loqL) = 1 
is a sharp turning point between propagating and evanes- 
cent waves. 

It might seem difficult to evaluate the summation Ij24(l 
analytically. But since we are solving exactly the same 
problem as the previous section, it should be clear that 
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[Cf. Eq. CU)] 
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E^ 2 = 




(A0) 1D (/ P /A), (26) 



as we have verified numerically in the cases loqL = 50, 
100 and 200. 

Having obtained the cut-off length scale of A, it is easy 
for us to offer a simple explanation to the extra sup- 
pression factor of (Ip/X) using Fourier optics. In each 
transverse direction, the spatial spectrum of refractive- 
index fluctuation is flat, up to the spatial frequency of 
1/Zp. Since only fluctuations with spatial frequencies be- 
low 1/A propagates, we have a suppression of Ip/X in 
"power" , or \Jlp/\ in linear amplitude. With two trans- 
verse directions, we then expect the suppression factor of 
ip/A. 

Discussions and Summary. In this paper, we have cal- 
culated fluctuations on the phase front of an extra galac- 
tic point source, caused by Planck-Scale fluctuations in 
refractive index — a toy model for space-time foam. If 
diffraction of light were ignored, or if we assumed a space- 
time with one time dimension plus one spatial dimension, 
our toy model would give comparable results toprevious 
esimates on the random- walk model Q, 0, H3, Q| ■ How- 
ever, the diffraction of light forces us to average space- 
time fluctuations over all different possible optical paths 
that extend to all three spatial dimensions. This averag- 
ing filters out all fluctuations with transverse scales finer 
than the wavelength. In our model, this causes an extra 
suppression factor of Ip/X, with [Cf. Eq. JTJ] 



< ^/IpL/X 2 (Zp/A) 



This suppression is at least by 19 orders of magnitude 
for astronomical observations, if the wavelength of A = 
10~ 16 m were to be used. Numerically, we have 

A0 ~ 10~ 2 VV Gpc (5 x 10- 7 m/A) 2 , (28) 



which unfortunately leaves us with no hope of detecting 
this random-walk model with any conceivable astronom- 
ical observations. 



Even though in this paper we have only studied a 
static, random-walk-typc model in which space-time fluc- 
tuations are local to each spatial region with size Zp, 
we believe the diffraction of light is a general enough 
phenomenon to affect significantly the use of photons to 
probe any other model of quantum gravity. 



This work is supported by the Alexander von 
Humboldt Foundation's Sofja Kovalevskaja Programme 
(funded by the German Federal Ministry of Education 
and Research). 



Y.J. Ng, W.A. Christiansen, and H. van Dam, Astrophys. 
J. 591, L87 (2003). 

R. Lieu and L.W. Hillman, Astrophys. J. 585 L77 (2003). 
R. Ragazzoni, M. Turatto and W. Gaessler, Astrophysis. 
J. 587 LI (2003). 

W.A. Christiansen, Y.J. Ng, and H. van Dam, Phys. Rev. 
Lett. 96, 051301 (2006). 

S. Lloyd and Y.J. Ng, Sci. Am. 291, 52 (2004). 
V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 

1330 (2004). 

Michael Maziashvili, hep-ph/0605146 
J.D. Jackson, Electrodynamics, 3rd ed., John Wiley & 
Sons, New York (1999), Sec. 10.5. 

(27) 

[9] In order to obtain Eq. we break the right-hand side of 
Eq. IllCit into an integral over 5-functions: J dr' F(r)5(r — 
r'). For each S(r — r'), its contribution must be of the 
form Aje(uJor) for r < r' (due to regularity at origin) and 
Bhi (ujor) for r > r' (due to outgoing condition at infin- 
ity). Then A and B can be solved by applying continuity 
in ip, and junction condition in dtp /dr. 
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